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SUMMARY 


An investigation v?as made of various available analytical 
methods for predicting the aerodynamic performance of film 
cooled turbine blades. The methods were cataloged according to 
their basic approaches to the film cooling problem. Each 
method is briefly described, and an attempt is made to improve 
upon one of the methods. 

The improvements to the method involved the development of 
a technique to calculate cascade body surface pressure distri- 
butions with the effects of coolant injection on these 
distributions included. This is accomplished by using an 
integral equation solution for the potential cascade flow. The 
solution involves using a minimization criteria on the singularity 
density distributions which in turn are used to calculate the 
potential flow. The minimization technique was originally 
developed for the isolated airfoil problem and is extended in 
this report to cascades. Comparisons of calculated cascade blade 
surface pressure distributions with other analytical solutions 
and experimental data are made , and good agreement is shown . 

The modified method, which includes the coolant and primary 
flow interaction effects, was used to predict 'the aerodynamic 
performance of a film cooled turbine cascade for which cold flow 
experimental data was available. The method is in good agreement 
with experimental results , but offers little improvement over the 
performance predictions of previous methods. However some insight 
is gained as to the effects of the coolant injection on the total 
flow field such as interactions with primary flow, interactions 
between coolant rows,, injection induced separations, coolant mass 
flow distribution in full film cooling, and the onset of sonic 
injection conditions. 

Finally the method is used to make an analysis of a possible 
film cooling configuration. This analysis demonstrates the capa- 
bilities of the method for use in the design of film cooled 
turbine blades. 
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INTRODUCTION 


Background 

With, the increased use of turbofan and turboprop engines, 
a need has developed for higher turbine inlet temperatures . 

These higher turbine temperatures produce higher energy levels 
in the overall thermodynamic cycle so that more efficient 
engine designs can be made. This additional energy is also 
required to power the fan and prop components of these engines . 
Film cooling is one of the better methods of increasing turbine 
inlet temperatures. The film cooling process consists of 
injecting a coolant into the. mainstream so that it forms a layer 
of cool fluid which protects the turbine components from the 
hot combustion gases. Although injection or film cooling is 
known to be an effective means of allowing increases in turbine 
inlet temperatures, little is known about the effect it has on 
turbine aerodynamic performance, and even less about the local 
fluid dynamics of the injection process. In the present work, 
an analytical method is developed to predict these perfoinnance 
effects and to model the local inviscid fluid dynamic of the 
coolant injection process. The analysis is then used to study' 
the effects of film cooling on' turbine aerodynamics . 

The complexity of the actual turbine film cooling problem 
prevents one from solving it directly. The combination of three 
dimensional, unsteady and rotational fluid flow effects plus 
the heat transfer betv/een the cooling fluid and blade surfaces 
cannot be included in a single solution. In order to produce a 
manageable model of film cooling, several simplifying assumptions 
must be made. The present study deals with steady inviscid 
irrotational adiabatic flow through a two-dimensional cascade 
with injection. The solution also includes *the effects of 
compressibility, viscous losses, and the flow characteristics of 
the coolant ports. Figure 1 shows the basic parameters of the 
problem. The analysis is developed for known inlet and internal 
blade flow conditions. The downstream flow conditions and 
injection velocities are determined during the course of the 
solution. 
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Literature Review 


Before analyzing the aerodynamic perfoianance of film cooled 
turbine cascades, one has to gain an appreciation of the various 
basic techniques for the aerodynamic evaluation of turbo- 
machinery blades. References [1] and [2] are two general works 
that- treat the aerodynamic design of turbomachines. Stewart's [3] 
method is the classic analysis used by most designers in 
evaluating the performance of cascades. This method can be used 
in a totally analytical study as it was originally intended, 
or to reduce experimental data as describedby Prust [4]. 

The number of film cooling studies was limited before the 
late 1960's. References [5] through [8] represent a cross- 
section of the work presented during that time. Also text 
books such as Schlichting [9] and Kays [10] have sections on the 
general problem of fluid injection into a main stream. From 
the late 1960 's to the present, a considerable amount of research 
has been conducted on the film cooling process . Much of this 
research was instigated by the need for high temperature turbine 
components. References [11] and [12] give comprehensive 
information on the turbine cooling problem. 

Few studies concentrated on the aerodynamic losses associated 
with injection processes. References -[13] through [16] represent 
reports from NASA Lewis Research Center's on-going experimental 
program. Hartzel [17] and Hiroki and Katsumata [18] investigated 
experimentally turbine cascade aerodynamic losses as well as 
heat transfer characteristics. An experimental study of an 
actual turbine is presented by Lokai and Kumirou [19] , but the 
predominant coolant discharge was radially out of the blades 
end plate and from the blades trailing edge. In -the analytical 
studies of Prust [20] and Hartzel [21] mass averaging mixing 
theories were used to resolve the film cooled cascade performance 
problem. Tabakoff and Hamed [22] approached the problem using 
boundary layer concepts . 

A number- of analytical studies investigated the effect of 
injection on the boundary layer of the primary flow. While 
references [23] through [28] were concerned with the laminar 
boundary layer problem. Herring [29] analyzed turbulent boundary 
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layers with injection. Experimental data to complement Herring's 
analysis can be found in reference [30] , Most of the boundairy 
layer studies listed above included comments on heat transfer. 
Goldstein [31] and Ekert [32] describe some of the analytical 
aspects of the heat transfer problem. Experimental work on 
the same sxibject can be found in references [33] through [41] . 

Since potential flow integral equation solutions are used 
in the analysis of this report, some of the available studies 
in the field will be discussed. A review of the various tech- 
niques for solving the integral equation problem is given by 
Hess [42] . Two methods of solution have evolved. One uses a 
vortex singularity distribution as described by Martensen [43] , 
and the other uses a combination of vortex and source-sink 
singularity distributions [44] . In this study the second 
method is employed, but examples of the use of Martensen method 
can.be found in references [45] and [46]. The source-sink and 
vortex combination method has been continually developed and 
refined. Hess [47] described a method for obtaining higher 
order solutions. Geising [48] and Geising and Smith [49] extended 
the method to cascades and hydrofoils, Bristow [50] presented 
a singularity density minimization technique which produced a 
higher accuracy first order method. A rapid solution method 
was demonstrated by Dilley [51] for possible use in interacting 
boxandary layer solutions . A boundary layer corrected potential 
solution for multi-component airfoils using vortex distributions 
is reported in reference [52] . 

The turbine blade shape with its blunt trailing edge poses 
a question as to how to apply the Kutta condition. Most methods, 
[42] and [50] , simply require that the velocity on the upper and 
lower surfaces be equal , or that the noannal velocity to the 
trailing edge bisector be zero. Gastelow [53] , [54] and 
Geller [55] offer other possible solutions to the blunt trailing 
edge problem. 
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REVIEW OF AVAILABLE ANALYSES 


A brief description of the various solution approaches 
studied is presented. This provides background for the final 
solution method used in the present work, 

Boundaary Layer Methods 

The boundary layer approach has been used in the past to 

give loss estimates for cascades without injection. Several 

solutions have been developed which use a pressure distribution 

as a driving function for the boundary equations [56] , [57] . 

However, this standard approach fails when applied to the 

film cooling problem. The Prandtl boundary layer theory 

assTmption that, the dimensionless velocity in the normal 

1/2 

direction to the outer flow remain of order (1/Re) * , is violated 

by the injection process. Solution of the boundary layer equations 
with injection has been obtained [24] and [25] , when injection 
velocities were kept very small. Even with the injection 
velocity maintained at the proper level, the boundaiy layer can 
be induced to separate by the injection. The numerical 
solution usually fails in such cases of. produces, erratic results , 
The boundary layer approach has been applied to film cooled 
cascades with limited success [22] .• In this solution, the 
injection velocities were made small by applying the injection 
over a larger area than existed in the real case. Furthermore, 
the solution only produced useful results for injection from 
the forward portions of the blade where the flow is accelerated 
causing a favorable pressure gradient. An at.tempt was made to 
extend this method by including the inviscid effect of the 
injection on the outer flow. This inviscid effect was noted 
experimentally by Bergeles , Gosmann, and Launder [37]. The 
extension allowed calculation of higher injection velocity flows, 
but useful solutions were still limited to the forward regions 
of the blade. 

A study of the available interacting boundary layer tech- 
niques for solving the small or closed separation flow problem 
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points out some of the complexities of the injection process. 

Such a study could be based on ideas described in references 
I58j and [59]. In' which case, the requirement of a small 
disturbance due to the inj.ection would limit the study's 
application in the real film cooling problem. Also, the' finite 
injection area introduces two discontinuities in the flow field 
at its beginning and end. These discontinuities must be taken 
into consideration in the problem formulation. Furthermore, 
if the flow field is subsonic, the effect of the main flow over 
the port on the injection velocity profile should be considered. 

In conclusion, it can be said that present boundary layer 
methods are not capable of handling the turbine film cooling 
problem with injection velocities used in actual applications. 
These methods should prove useful in the transpiration cooled 
turbine, where the cooling flow is distributed over the entire 
blade surface. This will keep injection velocities small, 
eliminate the discontinuities in the boundary conditions, and 
remove any injection induced separations. The solution will 
still not be straightforward, since the injection velocity is 
a fxinction of the pressure distribution and is not prescribed. 
Also, if interacting boundary layer techniques are to be used, 
solution branching must be resolved by providing a downstream 
boundary condition for the boiandary layer equations [58] . 

Mixing Methods 

Mixing methods use a simplified analysis of the film cooling 
problem. The methods follow roughly the concepts of Stewart's 
[3] downstream mixing model. They are based on conservation 
laws and provide gross information rather than details of the 
injection process. This information is sufficient, however, 
to predict the aerodynamic loss of the cascade. Four mixing 
methods are used at the present time to provide estimates of 
turbine blade performance. Since these methods are described 
in references [20], [21] and [60], only the portions used in 
the present analysis will be discussed here. 
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Constant static pressure miKing . The simplest of the 
methods is the constant static pressure mixing model of reference 
[21]'. The cooling jet row is modeled as a slot inclined to the 
surface. Assuming a perfect gas, the equations of conservation 
of mass energy and streamwise momentum, are applied to the 
system depicted in Figure 2, The outer boundary of the system 
is free to expand and contract so that the static pressure 
remains constant. The momentxim of the jet normal to the primary 
flov7 is not considered in the analysis. 

The governing equations are: 

Conservation of mass 



u 


1 + 5 


Conservation of energy 



1 

1+5 


(T ■ + 5T^ ) 

u c 


( 1 ) 


( 2 ) 


Conservation of streamwise momentum 


V 


m 


1 

1 +? 


(v^ + 5V^ cose^) 




Ideal gas laws 

H. = T. 

t P t 

(4) 

P = p RT 

Where the different pressures, temperatures and velocities are 
defined in Figure 2 , and 5 is the ratio of the coolant to primary 
mass flow, W^^/W^. Note that the specific heat, c^, has been 
assumed constant throughout the flow field. 

The change in total pressure, AP^, between the upstream and 
mixed condition is determined using isentropic relations , and 
the assumption of constant static pressure. 
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M 


u 


Y 

Y-1 


where M is the Mach number and y the ratio of specific heats. 
Again using the constant static pressure assumption, the Mach 
n\3inber of the upstream and coolant flows are calculated from the 
local pressure distribution on the outer boiindary at the 
injection point 

Y-1 
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( 6 ) 



2 - 

Y-1 



Y-1 
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C 


1] 


(7) 


These Mach nximbers are also used to determine the upstream and 
coolant velocities, and V^. 

The mixed flow Mach number may be calculated using the 
definition of critical velocity ratio. 


(^) = [ 

'‘V ^o '■ 

cr 




m 


2yR T. 


m 


Yt^*l 


( 8 ) 


and the relation between Mach nximber and critical velocity ratio. 


2 

m 


2 

Y+"l 




'V 


cr' 


2 

m 


1 


Y-1 

Y-J-1 


(^) 

cr 


2 

m 


(9) 


Combining equations (2) , (3) , (8) and (9) , the mixed flow Mach 
number is expressed in terms of the upstream primary and coolant 
flow conditions as follows ; 


2 

m 


(y-1) c^(1-?)(T + ?T, ) 

^ ^U ^C _ Y-1 ] 

- (v^ + COS0)^ ^ 


(10) 
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Using the above relations for Mach numbers and the pressure 
distribution over the turbine blade, the total pressure change 
across each injection port may be calculated. The total loss 
due to the injection is the sum of the individual total pressure 
changes across each coolant, port. This loss is added to the 
viscous, trailing edge injection and trailing edge blockage 
losses to arrive at an estimate of the cascade performance. 

Composite mixing method . As presented in reference [21] , 
the method consists of repeated application of static pressure 
mixing within a layer. In terms of computational efficiency 
and overall accuracy, it represents the best aerodynamic 
performance estimator of film cooled turbine .blades to date. 

The flow problem is depicted in Figure 3. The flow field 
is divided into three regions, two of which are mixing layers 
along the suction and pressure surfaces of the blade. The 
free boundaries of these surface layers correspond to the outer 
boundary in the static pressure mixing theory. The mass flow 
through each layer is the same at all times, so that the layer 
has to be considered to expand in the out plane direction to 
accommodate the injected flow. The percentage of the total 
flow considered to be within the layer is somewhat arbitrary, 
since the final performance results are only mildly dependent 
on this choice. Choosing the single layer mass flow based on 
twice the total amount of the expected injection flow in this 
layer has worked well for this investigator. The third part 
of the flow field is the main or primary flow, which foinns an 
inviscid core in the flow passage between the mixing layers. 

The static pressure is considered constant across the mixing 
layer, and is controlled by the inviscid core flow pressure 
distribution which is taken as that on the blade surfaces without 
injection. 

Various modified forms of the static pressure mixing have 
to be used depending on the injection location. Near the blade's 
leading edge most of the fluid is moving in the direction of 
the inlet flow. This direction is therefore adopted for stream- 
wise momentum rather than the direction tangent to the blade 
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surface. Along the blade surfaces the mixing theory is applied 
as it was previously presented. The contents of the two layers 
are mixed with the in viscid core flow at the trailing edge 
plane/ .and all three flows are assumed to be moving in the 
streamwise direction which is dictated by the trailing edge 
angle. 

The aerodynamic performance of the blade row is calculated 
by considering the total pressure changes in the mixing layers. 
Starting at the leading edge, the change in total pressure is 
calculated as each injection port is encountered. Unlike the 
one dimensional mixing model, however, the upstream total 
pressure at a port is. dependent on the previous total pressure 
changes in the layer due to upstream injection ports. At the 
trailing edge plane the layer and core flows are mixed using 
the constant static pressure method again. This provides a 
calculation of the change in total pressure through the blade 
passage due to the coolant injection. With this flow, a 
modified version of Stewart's method [3] which includes 
trailing edge injection, is used to complete the performance 
calculations . 

The modified Stewart's method was developed by applying 
the general analysis of reference [3] to the flow depicted in 
Figure 4. Trailing edge injection was included in the method 
by incorporating additional terms in the conservation 
equations. The final equations for the downstream critical 
velocity ratio are given below. 

Axial critical velocity ratio is 


V 

^cr 2 


yC 

Y+1 



2 


1 + 


T~1 

Y+1 


V 


'V 


cr 


( 11 ) 


where 
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Tangential critical velocity ratio is 


cr te 
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V - 
cr 2 


sina„[(l~d,-e^} (^) + A (^) (^) (^) ] 

2 11 te p T 


cr 


t te 


1 cr te 


,V 
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cr 1 '^t te cr te 


-) 


(13) 


In the above relations, is critical velocity ratio, 

is the area/unit span of the trailing edge injection slot, 5 
and 0 are the boundary layer displacement and momentum 
thicknesses, ct is the flow angle, and p is the fluid density. 

The subscripts 1, te, and 2 indicate flow properties at the 
trailing edge plane, trailing edge injection slot, and far 
downstream respectively. 

Other Methods 

Two methods for predicting aerodynamic losses were developed 
during the course of this analysis , but neither method produced 
acceptable performance estimates. They are presented, because 
some of their concepts are used in the analytical method of 
this report. 

The first method was an inviscid flow analysis that included 
inviscid injection effects. The cascade loss was calculated 
by perfoimiing a numerical integration of the flow properties 
downstream of the blade row. This method did not succeed because 
the inviscid flow alone did not contain enough of the physics of 
the loss generating mechanisms . The mixing and viscous proper i_ies 
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of the fluid proved to be significant contributors to the losses, 
so that the performance estimates were characteristically too 
high. 

The second method was suggested by Prof. Douglas E. Abbott 
of Purdue University, Lafayette, Indiana. He has had some 
tinpublished success with it in predicting aerodynamic perfor- 
mance for single airfoils. The method is the same as used by 
Geller [45, 55] for predicting cascade flows. The basic 
analysis simulates the trailing edge flow separation as an 
injection of fluid out of the blade surface. The injection 
produces a constant pressure slip line over the separation. 

The simulated separation flow is interacted with a boundary 
layer calculation to locate the separation point on the blade 
surface. Lift and drag forces are then found by integrating 
the inviscid flow pressure distribution and boundary layer skin 
friction along the blade surface, assuming the skin friction 
to be zero in the separated regions. 

When this technique was applied to the turbine cooling two 
problems arose. First, the coolant injection segments the 
boundairy layer, making any continuous boundary layer calculations ' 
unfeasible. A boundary layer solution reported by Beavers and 
Stratford [61] may be adapted to provide a local similarity 
skin friction estimator, but any injection boundary layer inter- 
action is ruled out if multiple ■ cooling flow injection is to 
be studied. However, the interaction process may be rendered 
unnecessary if the separation point is assumed to coincide with 
the injection port. This is supported by experimental evidence 
[37]. Second, difficulty was also encountered in determining 
the components of lift and form drag from the pressure distri- 
bution integration, since the primary flow is turned by the 
cascade, and selection of reference velocity on which to base the 
lift and drag vectors is not obvious. With the usual choice of 
potential flow solution onset velocity as a reference, cases 
were encountered where positive form drag was calculated. The 
method still has some merit as an alternative to mixing theory, 
but further study is needed. 
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IMPROVEMENTS TO COMPOSITE MIXING MODEL 


Of the various analytical methods studied, the composite 
mixing theory was found to give the best aerodynamic performance 
estimates for the widest range of injection locations and 
velocities. This method was therefore, selected as a basis 
for analysis of this report. 

Improvements to the composite mixing method can be obtained 
by modifying two basic assumptions of the mixing theory. These 
ass\imptions are that the in viscid core flow which drives the 
mixing theory is unaffected by injection, and that the effects of 
injection on the boundary layer flow are negligible. By using 
the inviscid flow solution developed here, the injection effects 
were included in the mixing method. The detailed structure of 
the boundary layer with injection was found to be too complex 
to lend itself to analysis. However, by using a method similar 
to Geller [45] , the effect of separation on the inviscid flow 
were modeled. 

The analytical method used in this investigation was developed 
by incorporating these improvements into the mixing theory. 

The present method, therefore, consists of the composite mixing 
theory driven by an inviscid flow which includes the effects of 
injection and injection induced boundary layer separation, 

POTENTIAL CASCADE FLOW PROBLEM FORMULATION 

’ USING DISTRIBUTED SINGULARITIES 

To incorporate the improvements to the composite mixing theory 
described in the previous section, a means was needed to calculate 
an inviscid cascade flow with injection effects included. The 
integral- equation or distributed singularities solution was 
chosen for this calculation. It was selected because the solu- 
tion can easily be adapted to calculate the injection effects 
on the inviscid flow. In addition the distributed singularities 
method has the potential to provide a very rapid solution since 
the cascade blade surface velocity distributions can be determined 
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without solving the complete flow field. References [42] through 
[49] and [65]/ provide detailed descriptions of the basic solu- 
tion derivation, therefore only a brief outline will be given 
in this report,. 


General Cascade Potential Flow Problem 

Referring to Figure 5, the cascade problem differs from the 

isolated airfoil problem in that the upstream and downstream 

velocities CV_ and V_) , are not the same. This is due- to the 
X h 

flow circulation induced by the cascade across the entire flow 
field. However, the potential flow problem for the two 
dimensional cascade does follow the same line of though as that 
for a single airfoil with two modifications. First, the 
contributions of each of the airfoils making up the cascade 
have to be considered. Secondly, in the cascade problem, there 
is no physical counter part to the invariant free stream velocity 
in the single airfoil problem, but the idea of an undisturbed 
velocity is used to define an onset velocity which is constant 
throughout the flow field. This velocity is determined during 
the course of the iterative solution of the present method. 

In the potential problem it is assumed that the flow is 
inviscid, steady, irrotational , and incompressible. With these 
assumptions the governing equations take the following form. 

Continuity is , 

V « V = 0 cm 

Conservation of momentum becomes simply Bernoulli’s equation, 
and is used to determine pressures in the flow field. 


P-P 

onset 
1 „2 
2 *^^onset 


= 1 - (: 


V 


C12) 


onset 


where c is the pressure coefficient. 

P 

Boundary conditions at the cascade inlet and exit are. 


V V, 


as 


R ^ 


upstream 


(13) 
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and V ^ Vg as R “ downstream (14) 

where R is the position vector from the cascade row. The other 
boundary condition is at the cascade body surfaces and specifies 
the normal velocity along these surfaces, F(t), 


V ' 



F(t) at R = 


(15) 


Solution of the Governing Equations 

The flow velocity , V, can be expressed as the summation of 
a constant onset velocity, ^onset' ^ variable disturbance 

velocity, v, which is not necessarily small 


V 


V 


onset 


+ V 


(16) 


Since the onset velocity is constant throughout the flow field, 
the momentum is satisfied if the disturbance velocity v, is 
described by the velocity potential function ij) as follows. 


V = - V(J> 


(17) 


By substitution into equation (16) , the velocity vector becomes 


V = V . “ V(j) 
onset 


(18) 


Using this definition in the continuity equation (11), Laplace's 
equation for the potential function is obtained 


V (|) = 0 


(19) 


There are a number of real harmonic functions which satisfy 
the above equation, but it is more convenient to work with 
complex variables and take advantage of their compact notation 
and special algebraic manipulations. 

The potential function (<j)) is the"real portion of the 
complex potential, C(z) , 

CCz) == (z) + i i(j (z) (20) 
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where z is a point in the complex plane 

z — X + iy (21) 

In the above equations, i denotes the imaginary number (/^) , 
and Tp is the stream function. The coordinates of z in the 
complex plane are x and y. An analytic function may now be 
selected for G(z) , so that the real part ((j>) will satisfy 
Laplace’s equation. 

From ideal flow theory, two singular functions, namely 
source-sink and vortex, are known to be analytic throughout the 
flow field except at their origins . If c =s a + bi is the 
singularity location and r is the distance between c and z, 
the point of interest in the flow, these functions will have 
the following forms. The source-sink velocity potential is 
given as 

4 Cz-c) = Re iCCz-c)] = s An r (22) 

and the vortex velocity potential is written as ; 

^ (z-c) = Re [C(z-c)3 = - g8 (23) 

where polar notation is used for the complex argument of the 
functions., so that 

z-c = r e^® (24) 

In the above equations, g and s are defined as the vortex and 
source-sink strengths. These potential functions may be written 
as a single function by defining a combined singularity strength. 

k = (s + ig) (25) 

The complex potential then becomes 

C(z-c) = k An(re^®) (26) 
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and the combined velocity potential is: 


(j)(z-c) = s An(r) “ g9 (27) 

Laplace's equation is linear, and the sum of analytic function 
is also analytic. Using these facts, a solution may be formed by 
superposition of several singularity functions located throughout 
the flow region. For the single airfoil, the singularity functions 
are placed along the body surface to produce the disturbance 
potential, but for a cascade the simple sink-source and vortex 
potentials do not account for total effect of all the airfoils in 
the cascade row. 

Referring to Figure 6 , the complex y axis is aligned with the 
blade row, and the real x axis is perpendicular to it. Since the 
cascade is infinite, corresponding singularities on the cascade 
bodies located at distances S parallel to the y axis will have 
the same strength or density, k. The complex potential at any 
point in the flow field due to a line of these singularities can 
be written in the following series 

CCz) = kEAn(z-c) + £n(Z“C+iS) + 5,n(z-c-iS) + An(z-c+2iS) 

+ ...] (28) 

Combining. terms v?ith the same multiples of iS the series becomes 

C(z) = k{£n(z-c) + £n [ (z-c) + £n [ (z-c) (2S) ^3 

+ . . . } (29) 

The components of the disturbance velocity Cv^rV^) are found by 
taking the derivative of the complex potential. With a rearrange- 
ment of terras the resulting series takes the following form; 

dC . _ ^ f S 2-iT Ca-c) - 1 

dz S ^tTCz-c) S ^2 2, ^2.2 

IT +7T (z-c) /S 


+ K 5 + .».] (30) 

4Tr^+r'^(z-c)^/S 
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which is the same as the defining series for the complex hyper' 
bolic cotangent. 

The complex disturbance velocity, w, can now be written 
as a single functipn,. 


V - 
X 


/Tf (Z”C) 1 
s~ g ) 


and the complex potential may be found by integration. 


C(z) = k An [sinh ] (32) 

This is the elementary potential for a cascade of singularities. 

A body can be represented by distributing the source-sink 
and vortex singularities along its surface. The simplest way 
to distribute the singularities is to replace the body with 
polyhedron and use a line singularity of constant density over 
each segment of the polyhedron. The basic geometiry of the problem 
is illustrated in Figure 6. To develop the line singularity 
formulation for the velocity function, one integrates the 
elementary velocity functions along the line segment on which 
they are distributed. A general segment lies along a coordinate 
line C oriented at an angle n to the x axis of the reference 
complex plane. The integration is between the bounds, of the 
line segment c^^ and C 2 » This discussion has concentrated on 
one body, but similar statements may be made for a cascade of 
bodies so that the integration takes place on the complex series 
derived previously in equation (31) . The complex flow disturbance 
velocity, W, for a cascade can be expressed as follows: 

^2 

W = / coth( ^'^f~^^ ) dc (33) 

* c ^ 

^1 

where 


c 

= ih + ? e ^ 

f 

(33a) 

dc 

= e^’^ dc 

f and 

(33b) 


= e dc 


(33c) 
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then 


W(z) 


kir 

S 


(34) 


f 

c- 


coth 




dc 


Integration gives the following form of the complex velocity; 


W(z) 


sinh [it (Z-C 2 ) /S] 
^sinh tir (2-c^) /S] ^ 


(35) 


The final form of the complex velocity for the two-dimensional 
problem requires another integration of the line singularity 
function in the out of plane direction to plus and minus 
infinity. This integration results in a factor of 2 and repre- 
sents a plane surface singularity distribution: 


W(z) 



sinh [ it (z~C 2 ) /S] 
^sinh [it (z-c^) /S] ^ 


(36) 



With the disturbance velocity given by equation (36) , which 
satisfies the governing equation, it is necessary to have a 
closer look at the far upstream and far downstream flows in 
order not to violate the cascade inlet and exit boundary condi- 
tions. At a far stream station, the cascade singularity 
distribution would appear as a row of point sink-source and 
vortex singularities. 

The disturbance ■ velocity at the far upstream and downstream 
boundary conditions is determined by taking the limit of equation 
(31) as 


^ (z-c) -5- t '•«> far downstream 

W 

•— (z-c) ^ - 00 far upstream 


(37) 


and by replacing kby K=E+ri in which S and r are the 
concentrated source-sink and vortex density of a cascade blade. 
After taking the limit, the upstream disturbance velocity is 
found to be 
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X 


= 0 


rir 

s 


(38) 


Similarly the downstream disturbance velocity is given as 
= 0 


V 


X- 


V 


+05 


+C0 


Tit 

s 


(39) 


Now from ideal fluid floW/ the circulation, r . , is given by 
the line integral of the tangential velocity component around 
the body and can be related to the vortex density as follows: 


cxr 


= .^ V • dt = 2irr 


(40) 


Using equation (40) in equations (38) and (39) , the following 
relationships are obtained. 


V 


V 


cir 


2S 


cir 




2S 


(41) 

(42) 


These disturbance velocity components, snd 


are shown 
respectively 

Using these diagrams , the following relations can be 
written; 


in the vector diagrams of Figure 5 , as V and V, 

^ ' up down 


a . 
xnc 


3 


(43) 


A0 = 


“l “ “e 


(44) 


where a. is the incidence angle of attack, a_ is the inlet flow 
me X 

angle, and 8 the stagger angle of the cascade. Also A0 is the 
flow turning angle and is the exit flow angle. 

Using equations ^(40) and (41) for and and the 

vector diagrams of Figure 5, the following relations are obtained 


= tan - ['^onsat 2^"" 

. cos a •• 

onset 


-1 


(45) 
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r . 
cir 


•] 


( 46 ) 


1 sina - 

, “1 r onset 2S 

a_ == tan I — ST 

E V L COSa 


onset 


where a is the angle of attack of V^^set negative as 

shown in Figure 5. 

In the solution being presented, inlet flow conditions are 

given. Equation (45) is then used iteratively with the velocity 

solution, to determine the values of V a, and r^. . Once 

onsei- cir 

a velocity solution which is consistent with the inlet conditions, 
is obtained, equations (45) , (42) and (43) are used to calculate 

a„, A6, and a. . Note that if both a_ and a„ are given the 
velocity solution parameters (a, and V^^set^ found 

without iteration. 


Resolving the Cascade Body Surface Boundary Condition 

The cascade body surface boundary condition is the final 
requirement to be satisfied by the solution. From this boundary 
condition comes the integral equation to be solved. Recalling 
equation (15) , the boiindary condition was written in general 
form as 


V* n I 


= F(t) 


at R = R] 


body 


Substituting equation (18) for V into the above equation the 
following relation is obtained 

’♦•"It = Vset • "It - 

Using the polyhedron approximation of the blade, control points 
are selected at the mid-points of each of the sides of the 
polyhedron as shown in Figure 7. The surface boundary condition 
is satisfied at each of these control points. The disturbance 
velocity at the control points is the total effect of all the 
surface singularities along the polyhedron sides . 

n 

v = v. -iv =5-V(|) = ^W(z) (48) 

^ ^ i=l 
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Substituting equation (36) for W gives; 


V 


n 


I 

m=l 


-2k e 


-in 


An 


3lnh[it(z-Ci^^l)/Sl 

'•sinhiir (z-c^) 7s] ■* 


(49) 


Substituting this into equation (47) , the following relation is 
obtained ; 


n 


I 

m=l 


-2k^(t) 


sinh[iT (z-c . , )/S] 

—in r Itl+1 

® ^Sinh[ir(z-cj7s] 


onset 


■n 


- F . (t) 
(50) 


Where j- is an index that identifies the control point at which 
the boundary condition is being evaluated. It should be noted 
that surface singularity of equation (36) was derived in a 
coordinate system aligned with its surface. When this relation 
is used in the integral equation (50) care must be taken to 
resolve the disturbance velocity components into the coordinate 
system of point in the flow field where the flow velocity is to 
be calculated. 

By writing relations similar to equation (50 ) , for each 
control point on the polyhedron-, a system of n equations will 
result. The solution to these equations gives the singularity 
density, k^ (t) , at the control points. 


kj (t) - Sj (t) + i (t) 


(51) 


The above derivation gives the discrete approximation of the 
integral equation. The integral equation form may be recovered 
may be recovered by using an infinite sided polyhedron to 
represent the blade: 

/ V k(t) f(t) • n(t) dt = V^nset * ^ 


where f(t) is the known singularity function and k(t) is to 
be determined. 

The value of the function F(t) is noonnally taken to be zero 
for the solid boundary case. If injection occurs, F(t) or F^ (t) 
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takes on the value of the normal velocity component of the 
fluid injected at the particular location on the blade surface. 
After the singularity density distrib hbion has been found, 
the velocity at any point in the flow field may be calculated 
using the following relation; 


+ iV 
X y 


V cosa + 

onset 


V ^ sina + v 
onset X 


+ rv 


(53) 


Where v^ and v^ are the components of the disturbance velocity 
V which can be determined from equation (49) . 


DISCUSSION OF SOLUTION TECHNIQUES FOR THE 
POTENTIAL FLOW PROBLEM 

Several solution techniques exist for developing the singu- 
larity distributions. In most of the methods, which use a 
combination of source-sink and vortex densities, the singularity 
distributions are considered to be independent of one another. 

The simplest techniques for finding the distributions assumes 
constant source-sink density over all the elements . The values 
of source-sink density for each element are found by solving the 
system of equations similar to equation (50) for each element. 

The vortex density is determined from the Kutta condition or a 
similar requirement on the flow leaving the trailing edge. 

This technique will mainly be used for comparisons and is there- 
fore, referred to as the base method. By using more complex, 
element geometries and singularity distributions, higher order 
accurate solutions than the base method may be formulated [47] , 

[51]. 

Bristow [50] developed a method for single multi-body airfoils. 
In his method the source-sink and vorticity distributions are 
found using a density minimization function, which makes them 
interdependent. Several improvements over the base method are 
accomplished by using this technique. Although both methods are 
first order in error resoltuion the minimized solution has 
better accuracy. The comparison problem of a single circular 
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cylinder without circulation was studied, and solution errors 
are compared at the peak velocity point 2 1/2 radii above the 
cylinder. As can be seen in Figure 8, the minimized method 
produces the same error level as. the base method with iess 
than half the number of elements being used to approximate the 
cylinder. Another advantage of the minimization solution is 
that it can easily handle very thin geometries. The base 
method cannot analyze these geometries without reducing element 
sdze in the thin region. Also the solution in high velocity 
gradient areas (i.e. airfoil leading edges) is improved by 
minimization. Finally, the base method is very sensitive to the 
relative size of adjacent elements. No element should be double 
the size of element next to it, or a spike in the solution 
curve will result at the element. The minimized method is less 
sensitive, but violation of the element size criteria will 
still produce erratic solutions. Because of its simplicity 
and improvements over the base method, Bristow’s technique 
was adapted to the cascade problem and used in the analytical 
study. 


SINGULARITY DENSITY MINIMIZATION SOLUTION FOR~CASCADES 

The minimization method uses an expanded form of equation 
(53) as a starting point. Two mutually perpendicular velocity 
components at any point i in the flow field are; 


V 

T. 

X 

^onset 

cos (6^-a) 

n 

+ I 

j=l 


n 

s. + I 

3 j=l 

*v. . 


(54) 

X 

“ ^onset 

sin (0 j^-a) 

n 

+ 1 
j=l 

B 

"ij 

n 

s f 1 
^ 1=1 

®v. . 
11 


(55) 


The angles to the reference coordinate system of and 

are 6^ and a, respectively. The source-sink singularity distri- 
bution,, Sj, and the vortex singularity distribution, g ^ , are to 
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be determined, however the general forms of these distributions 

must be specified before solution is initiated. Using these 

specified forms and the given geometry of the problem, the 

influence coefficients (A and B ) due to the source-sink singu- 

s s 

larities and the. influence coefficients (A.^ and B-^) due to 
the vortex singularities, may be determined. If equation (55) 
is written for the control points of the body, the system of 
equations that results in equivalent to the body surface boundary 
condition previously described by equation (50) with replacing 
F.(t). In matrix form, the system is expressed as: ^ 


[V ] 

^ N. , 
X nxl 


'-Vset 


+ [B_ ] 


nxl 


X] nxn 


[s ] 

^ nxl 


-h 


tB ] 

V. . 

xj nxn 




nxl 


(56) 


Solving for the source-sink strength distribution gives 


nxl xj nxn x nxl 


-1 

- [B ] IB ] [g.J 

X] nxn xj nxn nxl- 


(57) 


Defining the following new variables 


'V , = ’'onset 

nxl X 3 J nxn x nxl 


(58) 


-1 

nxn ij nxn x] nxn 


(59) 


and using these variables [equations (58) and (59) 3 / equation 
(57) can be written as follows; 




nxl 




nxl 




nxn 




nxl 


(60) 
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A density function which is to be minimized# can now be formed in 
which the only unknowns are the strengths of the vortex 
distribution. 


g^3 


n 


I 

k=l 





V, 



(61) 


where and A^ are the areas over which a particular singularity 
has influence. The source-sink strengths (Sj^) are eliminated 
with the use of equation (60) . 

In order to minimize the function G and also to evaluate 
the influence coefficients, the form of the singularity distri- 
bution must be known. The form of the source-sink singularity 
distribution, Sj^, is assumed to be simply a constant over each 
element. However,- since minimization involves- taking deriva- 
tives, the function, G, should be continuous and differentiable 
with respect to its independent variables, g^^. To assure this, 
the vortex strength distribution, gj^, is selected to be, linear 
over each element and continuous around the body. Referring to 
Figure 7, this implies that the final value of g on one element i 
is the initial value of g on the next element, i+1. Assuming 
the elements and body points are numbered as in Figure 7 (i.e. 
starting at the lower trailing edge and moving in a clockwise 
direction) , a discontinuity in the vortex strength distribution 
will occur for closed body at the trailing edge, unless the 
value of g^ is made equal to the value of This also helps 

to keep the solution bounded in the trailing edge region. 

With the forms of singularity strengths known, the influence 
coefficients and distribution areas may be calculated. Evaluation 
of the influence coefficients require considerable algebraic 
manipulation and lengthy equations , The derivation of coeffi- 
cients is therefore presented in an Appendix. The distribution 
areas, the A's in equation (61), are calculated as a fraction 
of the total perimeter of the body. If each element is long, 
the following relations can be written ; 

n 


A 




I 

k=l 


(62) 
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for k 1 



n 

Uy. + I S,^ 

k“X 

n 

iz, + Z )/2 J 5,, for k = 1 

^ ^ k=l ^ 


( 63 ) 


The final step of the solution is the incorporation of a 
requirement on how the flow is to leave the trailing edge of 
the body. For the sharp trailing edge airfoil shape, the Kutta 
condition, that the flow leave the trailing edge smoothly, is 
the obvious choice. However, a turbine blade has a round 
trailing edge, and no clear cut condition exists. In the 
solution being presented, the trailing, edge' geometry of the blunt 
ended airfoil is ignored. The problem is solved as an open 
body as depicted in Figure 7. This is justified by real flow 
considerations that the rounded trailing edge will be submerged 
in a viscous separation wake and therefore has little influence 
on the inviscid flow. For a trailing edge flow requirement, 
zero normal flow velocity is maintained on th e bisector of 
the trailing edge angle at a point downstream. This point is 
taken at a distance which is one- tenth the length of the smallest 
of the two elements adjacent to the trailing edge. This is the 
same criteria used by Bristow [50] , and is also compatible with 
the Kutta condition for sharp trailing edges. Physically, this 
condition may be construed as requiring a straight wake leaving 
the blade parallel to the trailing edge bisector. In order to 
■implement this condition, the solution must remain finite or 
bounded in the trailing edge region. To assure this, additional 
requirements are made that the values of the singularity strength 
on both sides of the trailing edge be equal. 

Mathematically the trailing edge conditions act as constraints 
on the singularity density function G. The zero nomal flow 
velocity requirement is written as i 


V, 


N, 


= -V 


TE 


onset 


sin (6^g-a) 


n 

I 

j=l 


B 


’TE- 


n 

I 


B 


V„ 


j=l TEj 


(64) 
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Siabstituting equation (60) for , equation (64) becomes 


-''onset 


n 


gj^) 


n 

+ I B 


j = l '^TE j ^ - k=l 

= 0 


j=l ^TEj 3 


( 65 ) 


or in matrix notation, it takes the following form, 

Vonset sin(0 -«) = [B ] [g ] 

onset TE s^j 3 


[B^ ] ([C.,] [g.] + [B ] Eg 3 > (66) 

“TEj Ixn nxn nxl TEj nxn nxl 


The trailing edge requirement for the source-sink singularity 
strengths is ; 

= s (67) 

1 n 

Using equation (60), equation (67) becomes, 


6, - B 
1 n 


^^nk “ .^Ik^ 

Ixn nxl 


(68) 


The same requirement when applied to the vortex singularity 
strengths is written as 

gi = gj^ (69) 

and was already enforced by the selection of a piecewise continuous 
vortex distribution. 

The problem can now be set up as the minimization of the 
density function, G, subject to two constraints. The method 
of Lagrangian multipliers is employed. 

Q = G + X^(eqn. 66) t ^2 

where and ^2 Lagrangian multipliers . Expanding 

equation (70) , the following is obtained. 
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Q = {[A^ ] UEB.] . + [C.,1 Eg,] )^3 

“i Ixn nxl nxn nxl nxl 


+ EA„ ] Eg?] } + a.EEb^ ] EB.] 

1 Ixn nxl ^TEj Ixn ^ nxl 


+ EB 3 EC 3 Eg^3 + Eb 3 Jg.] 

^TEj Ixn nxn nxl TEj Ixn nxl 

■ ''onset 

-C6i-S^)) (71) 

In functional form, equation (71) may be written as 
Q = QCg^A < 32 ' ^ 3 ' ••• ^n' ^1' ^2^ 

To find the minimum value of Q, derivatives are taken with respect 
to each of the independent variables and set equal to zero. 


= 0 


The above relations form a system of linear equations, and are 
solved for the value of the independent variables (g^ ... g^, X^, X 2 ) 
This gives the minimum. Q and the corresponding constrained 
singularity density distributions. 
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where D is a square coefficient matrix with dimension n+2 and 
E is column constant matrix with dimension n+2. Details of .the 
coefficient and constant matrices (D and E) are worked out in 
the -Appendix-, 

Analysis of equation (74) brings out some important facets 
of the formulation. The coefficient matrix, D, is dependent 
totally on the geometry of the problem. The constant matrix, E, 
is a function of the given values of the injection velocity 
over the surface and the angle of attack of the onset flow, a. 
Once the coefficient matrix has been inverted, an inviscid flow 
solution algorithm has essentially been created for the body 
geometry under study. Solutions for varying angles of attack 
and injection configurations become simply a matter of matrix 
multiplication, ■ 

The solution to equation (74) gives the vortex strength 
distribution. This is substituted in equation (60) which yields 
the source-sink strength distribution. These distributions 
are substituted into equation (54) using the control points on 
the body surface to give the tangential velocity on- the body. 

The tangential velocity is integrated around the body giving 
the circulation. The circulation is used in equation (45) so 
that the far upstream boundaory condition for the cascade problem 
is satisfied. The specified normal velocities and the calculated 
tangential velocities are used in the momentum equation formu- 
lation of the pressure coefficient to give the pressure distri- 
bution over the body. Finally, equations (54) , (55) and the 
pressure coefficient may be used to determine the velocity and 
pressure at off body points in the flow field. 

COMPRESSIBILITY^ CORBECTION 

The integral equation solution gives only the potential or 
ideal flow through the cascade. In a turbine blade row, local 
Mach numbers may approach unity because of the high camber of 
turbine airfoils , the reduction of flow area in the blade 
passages, and the favorable pressure gradient across the row. 
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Compressible effects become significant as a Mach nimber of 0,5 
is exceeded, so that a correction to the turbine potential 
flow solution is needed. 

Several techniques of varying levels of sophistication and 
computational effort are available to provide these corrections. 
The tangent gas or Kaimian-Tsien compressibility correction is 
one of the simpler methods, and was used in the solution being 
presented. The method is well known, and any text on gas 
dynamics contains a discussion on it (i.e., see Shapiro [60]). 
The tangent gas correction is presented below as it was used. 

For the method a correction factor is formulated 


X 


M 


2 

ref 


(!+' /l-M 


2 

ref 


) 


(75.) 


M ^ is the Mach number of the flow which has the same thermo- 
ref 

dynamic properties at the point of tangency between the ideal 
equation of stats and an approximate linear equation of state. 
For this solution, 'was the Mach number of the potential 

flow onset velocity. The compressible velocity is found from 
the following relation ; 

(| ) a-x) 

,V ^ ref inc 

V >2 

ref comp 1 - A (== ) 

'^rsf inc 


where the subscripts comp and inc indicate compressible and 
incompressible. The incompressible velocity ratios are obtained 
from the potential flow solution. Critical velocity ratios are 
then foimd using the compressible velocity of equation (76) . 

The remaining flow properties can be determined by using isen- 
tropic relations. The compressible pressure coefficient is found 
to be; 



Y 

P-P-r , 1 ? -57 7*^1 

[{1+ M^d-(^) )} -1] 

i p Vj yMj ^ ^ 


(77) 
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The ratio of specific heats (y) of the actual flow is used in the 
relation/ although y=-l would be in accordance with the Ka3onan- 
Tsein correction. The actual flow specific heat ratio, however, 
gave results which were in better agreement with experimental 
findings . 


INJECTION PLOW 


A method of calculating the injection velocity was needed. 

The problem was modelled as a total pressure reservoir discharging 
into a known static pressure through an orifice. The flow 
parameters are the same aS' those depicted in Figure 2 for the 
mixing theories. Using isentropic relations, the- injection- 
flow Mach ntimber is : 



2 

Y-1 
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(78) 


The injection Mach number, M , is used to calculate the injection 
critical velocity ratio as follows: 
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where 


cr 


2y 

Y+1 


R T. 


(80) 


and R is the gas constant of the coolant. 

Using the above relations the ideal coolant velocity, V , is 
determined. In order to obtain real flow velocities, experi- 
mentally deteimiined discharge coefficients, Cj^, for the injection 
row were used. 
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Areas equivalent to the total port areas of the injection 
rows are calculated so that injection mass flow estimates can 
be made using the two-dimensional model prob Jem. This is done 
by determining a slot width which gives' the same injection 
area per unit span as the row of ports. This slot width is then 
used in the analytical solution; 


SEPARATION MODEL 

■ The injection induced separation model, is schematically 
represented in Figure 9 . The model is similar to that used by 
Geller [45] and [55] , to approximate trailing edge separations 
in cascades. 

The separation region is considered to lie under a constant 
pressure slip line. This line starts at the injection port and 
continues downstream at the injection induced pressure level. 

The constant pressure condition is enforced by using pseudo- 
injections downstream of the coolant port. The potential flow 
pressure coefficient along the boundary of the separation is 


= 1 - (v; 
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N, 


+ v: 


(82) 
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Where V- 


N 


V, 


sep 


is the actual or pseudo-injection -velocity , and 


T 


sep 


is the tangential velocity along the surface which is 
primarily affected by the geometry of the flow field. The 
inviscid flow field velocity, is tangential to the body 

and separation regions at all times. The inviscid pressure 
coefficient can therefore be written as: 


= 1 - Vl 


xnv 


inv 


(83) 


Along 'the slip line of the separation region, the inviscid and 

separation pressures are equal, so tha-fc th e pseudo-in jection, 

V„ , can be detemined as : 
sep 


V' 
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sep 


V? - 
xnv T 


sep 


(84) 
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The pseudo-injection velocity, Vj^, will be real only if the 
injection induced pressure is less than the pressure due to 
the flow field geometry alone. Physically this implies that 
pressure due to the biade- passage geometry must be greater 
than the injection induced pressure in order to keep the 
separation region contained. If it is not, the separation 
region will expand into the inviscid stream. Therefore, a 
separation region cannot continually extend downstream in a 
locally favorable pressure gradient. 

Incorporation of the separation model requires alterations 
in the potential flow solution. First the model, which adds 
an iterative loop to the solution procedure, acts as- an 
additional constraint. This separation constraint is usually 
compatible with the basic solution except when the separation 
extends to the trailing edge. In that case an inconsistency 
may arise between the model and the straight trailing edge wake 
constraint of the potential solution. When this occurs, the 
solution converges to a solution which honors the separation 
model constraint. Second, if the pseudo-injection velocities 
are chosen to be directed outward from the cascade body surface, 
the mass flow build-up results in increased flow velocities. 
Since the pressure coefficient is dependent on the magnitude 
of velocity and not its direction, the pseudo-injection mass 
flow build-up problem was partially eliminated by alternating 
the direction of pseudo-injections into and out of the blade. 

The alternating pseudo-injections did not change the shape of 
the inviscid velocity distribution curve from the case when 
the pseudo-injections were all outward, but it did reduce the 
inviscid velocity levels. It was noticed that a solution 
instability was encountered sometimes , if the separation region 
extended to the next downstream injection port. Logic had to 
be added to the computer program to prevent this by insuring 
that the pseudo- injection just upstream of the second port was 
directed into the blade. 
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MASS PLOW CORRECTIONS 


Two mass flow corrections are required in the analytical 
solution. First in the real turbine cooling problem, the siim 
of the primary and injection mass flow will remain constant, 
since the injection flow is taken from the main flow in the 
compressor. To include this in the analysis , the calculated 
incompressible injection mass flow of the solution is used to 
correct the inlet flow condition. However, the injection 
velocity calculation and inlet flow correction lag the cascade 
blade pressure distribution calculation s o that iteration is 
necessary in order to produce consistent solution. The second 
mass flow correction is needed to alleviate discrepancies between 
the mass flow calculated using injection parameters and that 
indicated by an integration of the analytical flow conditions 
at a downstream station. These discrepancies arose because 
large approximating blade elements , which did not properly model 
the injection row, were used in the potential flow solution to 
reduce computational effort. The use of the large elements 
resulted in a higher mass flow and a correspondingly higher 
average velocity out of the blade row than should exist in 
actuality. However, the real injection row flow parameters 
were used to calculate the actual injection mass flows so that 
the potential flow solution mass flow, and average velocity 
could be corrected to their proper values. 

COMPARISON OF PRESENT SOLUTION WITH OTHER ANALYTICAL 
SOLUTIONS AND EXPERIMENTAL RESULTS 

The ability of the analytical solution to calculate cascade 
flows was investigated by performing a series of comparisons 
with other analytical solutions and experimental results. These 
comparisons were carried out as the solution was developed, 
so that weaknesses in the present solution could be isolated. 

First the potential or ideal cascade flow calculations 
were verified by using Gastelow’s exact solution [62] which is 
based on the use of conformal transformations. Figure 10 shows 
the; blade shape analyzed and approximating polyhedron, which 



was used in the analysis.- The vertices of the polyhedron were 
generated by patching two spline curve fits of the blade pressure 
and suction surfaces together with the ci-rcul-ar segment forming 
the leading edge. The pressure coefficients obtained using 
the two methods are shown in Figure 11. It can be observed 
that the present solution is in good agreement with the exact 
solution everywhere, except at three locations. Discrepancies 
did occur near the trailing edge and just aft of the leading 
edge on both the suction and pressure sides of the blades. The 
major contributor to these errors was found to be the curve 
fit approximating polyhedron. Referring back to Figure 10, 
it can be seen that the leading edge errors occur in the transi- 
tion region from small to large elements where some elements 
are nearly double the size of their neighbors . ' This violates 
the element sizing requirements of -the analytical method and 
can therefore, explain the errors in that region. In the 
trailing edge region, errors can be attributed to the failure 
of the blade cuorve fit approximation to capture the cusped 
trailing edge of the actual airfoil. 

Next, results from the present solution were compared to 
the experimental low speed cascade data of reference [63] . The 
cascade airfoil had a blunt leading edge, a symmetric thickness 
distribution with a maximum thickness ratio of 0.1, and a cusped 
trailing edge. Lift curves are shown in Figure 12 for two 
cascade solidities. Analytical and experimental surface pressure 
distributions are plotted in Figure 13. From these figures it 
can be concluded -that the agreement between experimental and 
analytical results is good. The pressure distribution obtained 
using the base singularity method is also shown in Figure 13. 

It is clear that the present minimization method gives better 
results than the base me-thod in the leading edge high velocity 
gradient region and in -the thin trailing edge area of the cascade 
airJoil. 

The final test of the present method was a cascade of typical 
turbine blades shapes as shown in .Figure 14. Using this blade 
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blade shape , the first comparison with experimental data of 
reference [64] is shown in Figure 15. The comparison shows 
that the present method does a fair to good job of predicting 
the cascade flow with the largest error occurring at the pressure 
side trailing edge. This error indicates that the trailing 
edge geometry needs to be better represented and/or that a 
different constraint on the trailing flow should be used. In 
order to explore the validity limit of the compressible flow 
correction, a comparison was made with reference [66] . The 
turbine blade shape was the same as in the previous comparison, 
but the cascade stagger and spacing are slightly different. 

Figure 16 shows that, as expected, the correction became less 
effective as Mach numbers increased. It was determined that 
the blade surface velocity should not exceed 85% of the 
critical velocity if good results are to be achieved. Analytical 
results were also obtained for the same cascade with 13,4 percent 
coolant mass flow injection and the' coolant total pressure 
equal to the primary flow. The change in the blade pressure 
distribution due to the injection is sho\i?n in Figure 17. 

Although no experimental data is available for direct comparison, 
the flow near the injection rows was found to be similar to 
the qualitative description of reference [37] . 

With confidence in cascade flow solution established, the 
validity of the total aerodynamic performance calculation pro- 
cedure was studied. The procedure was verified by comparing the 
analysis wi-fch the data from the experimental test program of 
reference [16] . Figures 17 and 18 show that the analytical 
solution is in good agreement with the experimental data. The 
solution appears to give better results for injection from the 
pressure side than from the suction side. The largest discrepancies 
between experimental and analytical results are observed for 
cases involving injection from either side of the blade near the 
trailing edge. As shown in Figure 20, the results obtained in 
using the present method are also compared to the composite 
mixing method of reference [21] , which does not consider injection 
and separation effects on the invisid stream. Figure 20 indicates 
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that the two methods have about the same error level when compared 
to experimental data [16] . Multiple injection computational 
results are compared with experimental data and the composite 
mixing model. In Figure- 21. This comparison shows the present 
method to be better than the composite mixing method at the lower 
levels of injection mass flow and decidedly worse for injection 
mass flows above 20 percent. These multiple injection results 
are primarily of academic interest only since the maximum injection 
mass flow in a real turbine is of the order of 10 percent. 

After the completion of this comparison study, a statement 
on the capabilities of the present method may be made. First 
the present method does a good job of predicting potential and 
compressible cascade flows with and without injection. Second, 
when used to estimate film cooled cascade aerodynamic losses, the 
present solution works well as long as high injection rates, 
sonic injection, and injection induced total separation (i.e. 
past the blade trailing edge) are avoided. Last, although the 
error levels of the present method are equal to those of the 
composite mixing model, the present analysis should predict 
injection mass flow levels better. Occurrences of undesirable 
flow phenomena such as injection induced separation, sonic 
injection, and injection-primary flow interaction which cause 
large cascade losses can also be predicted using the present 
analysis. Using this information, film cooled turbine blade 
designs may be improved. 


DESIGN STUDY USING THE PRESENT METHOD 

The final version of the analytic "solution can be used as 
a design tool. To demonstrate this and also to gain further 
insight into the aerodynamics of the film cooling process, a 
design analysis is made. Since a considerable amo\mt of data was 
already available for the turbine of reference [16] , a redesign 
of its cooling ports configuration was the problem considered. 

Cooling Configuration 

Port locations are determined by considering the surface 
pressure distribution, manufacturing capabilities, and heat transfer 
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characteristics. There are several reasons to avoid leading 
edge injection, although experimental and analytical results 
have shown reasonable aerodynamic performance. First leading 
edge injection is highly disruptive to the flow field, inducing 
separation and tending to extend the stagnation flow region 
over a larger area of the leading edge. This can be partially 
attributed to the use of perpendicular injection, which is due 
to the inability to construct ports with small injection angles 
in the leading edge region. Leading edge injection also requires 
a- high chamber pressure to produce outward injection velocities. 
Furthermore unless the port row .is positioned such that cooling 
flow is divided in tow directions , some portion of the leading 
edge will not be protected by the cooling film. The cooling of 
the blade leading edge can be- achieved using other methods such 
as jet impingement on the inner surface. 

Injection rows should not be located in areas where the 
local pressures are unfavorable, since separation is likely to 
result. Referring to Figure 17, it can be seen that locations 
between 5-10 percent of the pressure surface and 35-45 percent 
of the suction surface would be poor choices for port locations. 
Injection from the aft portion of the blade suction side should 
also be limited, since high injection velocities and performance 
losses can result due to the low pressures in that region. 

Two final design criteria are discussed in the following. 
First, injection rows should occur in pairs so that the ports may 
be staggered providing better film coverage. Second in the thin 
region of the trailing edge, manufacturing of the injection ports 
is difficult, and the region is also structurally weak. Pin 
cooling can be used in the region, and the cooling fluid injected 
out of the trailing edge. 

Using the above information, a turbine blade design was made 
and is depicted in Figure 22. Impingement cooling is used at 
the leading edge. Injection occurs between 15-30 percent and 
55-70 percent of the surface distance on both sides of the blade, 
using pairs of injection rows. These rows are modeled as slots 
which have widths equal to 0.00735 units of blade chord. Film 
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cooling injections are all vectored at the same angle to the 
local contour. A transpiration insert is used for the 55-70 
percent injection on the suction side. This insert has a low 
discharge coefficient which- keeps the' injection velocities down, 
and a large injection area which allows good film coverage. 
Transpiration injection occurs over the total area of the insert 
and is at 90 degrees. The trailing edge is cooled using the 
pin cooling and trailing edge injection combination. The trailing 
edge injection is from a slot 0.0055 units of chord wide and is 
directed along the trailing edge bisector. The film cooling, 
transpiration and trailing edge injection discharge coefficients 
are 0.7, 0.05, and 0.7 respectively. 

Parameter Variation 

The cooling blade design was analyzed for the effects of 
variation in injection mass flow, angle, and total temperature. 

The results are displayed in terms of primary air efficiency, 
which is equal to the ratio of actual kinetic energy output of 
the blade row to the ideal kinetic energy of primary flow alone. 

The thermodynamic efficiency, which represents the ratio of 
actual kinetic energy of the blade row to the summed ideal 
kinetic energies of the primary and injection flows was also 
calculated and displayed. 

Injection mass flow was varied by specifying different 
internal total pressures. The film cooling injection angle and 
coolant total temperature ratio were held constant at 30 degrees 
and 1.0 respectively. Figure 23 shows the results for this case. 
Positive or outward injection from all ■ coolant rows was not 
achieved until a coolant total pressure ratio of 0.99 was 
reached. Primary air efficiency shows a steady increase which 
seems to be proportional to the injection mass flow. Thermodynamic 
efficiency shows a slight decrease until the injection mass flow 
reaches 10 percent. It then starts to decrease at a greater rate. 
These results are a consequence of the formulation of the 
performance parameters . The primasry air efficiency contains a 
multiplier in the form of a ratio of total mass flow to primary 
mass flow. As injection increases, total mass flow increases and 
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primary mass flow decreases which causes the multiplying ratio 
to become larger. An- increase in the injection mass flow 
results in a decrease in the thermodynamic efficiency because 
the denominator, contains the ideal kinetic energy which 
increases with injection velocity, while the numerator remains 
relatively constant. 

Variations of the two efficiencies with injection angle 
is shown in Figure 24. The injection total pressure and tempera- 
ture ratios were set at 0.99 and 1.00. Both efficiencies decrease 
for increasing injection angle. This is due to the loss of the 
contribution of the injection flow to the streamwise momentum. 

With pressure ratio set at 0.99, and injection angle at 30 
degrees the effects of variations in the coolant total temperature 
were studied and the results are shown in Figure 24. Since heat 
transfer is not considered in the analysis, the temperature 
difference between the coolant and primary flow mainly affects 
the injection mass flow rats. As the temperature ratio decreases 
the injection velocity and mass flow decreases, so that the 
dependence of the performance parameters, on injection mass flow 
may again be used to explain these results. 

Finally, the blade pressure distribution for the design 
problem v^ith and without injection are displayed in Figure 26. 

The injection pressure distribution differs only slightly from the 
no injection case, and separation was not indicated anywhere 
along the blade surface. These improvements in the injection 
pressure distribution from that of Figure 17 are due to the 
selection of injection row locations and the lower injection' mass 
flows of the design problem. 

Several observations may be made on the overall results of the 
study. First, although paired injection rows provide better film 
coverage, the two injection rows were foxind to interact with the 
main stream in such a way as to produce injection velocities in 
the downstream row which are twice those in the first row. Second 
a reduction in the coolant total temperature ratio from 1.00 to 
0.50 resulted in reducing the injection mass flow from 8.5 to 6.6 
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percent for the coolant total pressure and injection angle under 
consideration. This fact has to be taken into consideration if 
the cold flow analysis is to be applied to a real problem. 

Third, if it can be avoided, injection should not occur before 
or at the throat of the cascade passage , since sonic velocities 
are likely to result. Last, for a given coolant flow condition, 
larger injection angles are more disruptive to the inviscid 
core flow. 

CONCLUSIONS 

The major findings of this analytical study are concerned with 
the aerodynamic performan'ce of film cooled turbine cascades. 

Several film cooling configuration parameters were found which 
significantly effect cascade performance. 

Injection from the suction side of the blade and particu- 
larly the aft portion of the suction surface towards the trailing 
edge produces the largest losses or worst perfoinnance. These 
losses are sometimes attributed to boundary layer separation, 
but the analytical injection induced separation model of this 
report shows the flow in the region to be fairly stable. Two 
possible causes for the suction side injection were analytically 
observed during the course of the study. First, high primary flow 
velocities are encountered along most of the suction surface , which 
increases the difference between injection and primary flow 
velocities. The larger energies consumed in accelerating the 
injection flow through this difference create the loss. Second, 
analytical calculations showed that the flow angle near the 
trailing edge plane was decreased by suction side injection. 
However, the calculated blade circulation varied only moderately, 
and the average downstream turning angle of the cascade was changed 
only slightly. Thus losses could be attributed to the energy 
used in returning the locally disturbed flow to the downstream 
angle dictated by the circulation. 

For a given injection mass flow, pressure surface injection 
calculations generally gave performance numbers which indicated 
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lower losses than those for suction side injection at the same 
mass flow. This appears to be a consequence of the slower 
surface velocities which decreases the di Serence between 
injection and surface velocities and a continually favorable 
.pressure gradient, which keeps boundaiy and film cooling layers 
thin. The analytical model did indicate that an injection 
induced separation bubble can occur on the pressure side if the 
injection row is located just aft of the blade leading edge 
where a pressure dip and recompression cause an unstable flow 
condition. Injection in this type of flow region should be 
avoided since a hot spot on the blade surface is likely 
developed at the separation. 

As disruptive to the flow field as leading edge injection 
appeared in the analytical solutions, it produces surprisingly 
low losses. Reasons for not using leading edge injection are 
the development of large regions of stagnation flow on leading 
edge surface between injection rows which were indicated by the 
analytical model and the high internal blade pressures needed to 
produce a given injection velocity. Although trailing edge 
injection is of no use for film cooling, it was the one 
injection location found, which improves cascade performance by 
filling the viscous wake with higher energy injection fluid. 
Trailing edge injection is therefore the best means of internal 
cooling fluid discharge out of the blade, 

When compared to location of injection, variation in 
injection angle appeared to have the second most significant 
effect on cascade losses. The closer the injection is to being 
perpendicular to the blade surface the larger the losses and the 
more likely separation is to occur. This is substantiated by the 
analytical solution of this report and the experimental work of 
others. In the transpiration cooling method, where perpendicular 
injection is employed, the injection velocity must be kept very 
low if the total injection mass flow, and cascade performance 
are to remain reasonable. 

The least significant of the injection configuration para- 
meters in the study was the injection row grouping. Pairing of 
injection rows so that ports may be staggered to provide better 
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film coverage did not excessively alter cascade performance. 

The combination of injection rows from different surface loca- 
tions was found to produce additive performance characteristics. 
One important finding of the analytical study which had not been 
noted experimentally is the effect of mutual interaction 
between injection rows, which are closely grouped together. 

At the downstream row, a low pressure region is created by the 
upstream injection so that its velocity and mass flow can be 
several times the upstream row values in the low injection mass 
flow case. Alternately the stronger injection of the downstream 
row creates a higher inviscid pressure over the upstream row, 
reducing its injection velocity. The worst interaction effects 
resulted at the low injection mass flow rates, since the down- 
stream injection can dominate the upstream row and induce such 
a high pressure that fluid is forced into the blade through the 
upstream row. 

In addition to effects of the coolant injection on cascade 
performance, another finding of the analytical study was the 
interaction effects of the coolant injection on the primary flow. 
For the typical case of unseparated flow over a single row, 
the main flow acts as though the injection flow is a biomp on 
the body surface. The flow decelerates as it approaches the 
column of injection fluid and accelerates as it moves over the 
injection row. The peak velocity is located just past the 
injection row. The larger the injection velocity the stronger 
this interaction. The analysis showed that separation does not 
generally occur with the injection mass flows used in turbine 
cooling and can usually be avoided by a good choice of injection 
sites. But if the flow does separate , the analytical solution 
indicates that the geometry of the blade passage or a downstream 
injection will normally terminate the separation region before 
it reaches the trailing edge. It was also shown that grouping 
of injection rows close together almost always results in 
separation between the rows. The final analytical observation on 
interaction effects is that primary flow reacts strongly to 
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injection located just before or at- the blade passage throat. 

The injection at the throat often forces the main flow into 
trans- and supersonic flow regimes. 

As one reviews the literature on film cooled turbines, four 
aerodynamic perfprmance parameters are found to be used by most 
investigators. After using all these parameters in this study, 
a ranking of their usefulness may be made. Of the four, the 
thermodynamic efficiency is the best indicator of injection 
cascade performance. Primary air efficiency's chief merit is that 
it clearly indicates the effect of -hh e coolant mass flow on 
performance. However, it often takes on values greater than one 
because of the ratio of the injection plus primary mass flow 
to primary mass flow used in its formulation. This may give 
the designer a false sense that he is getting more performance 
out of the injection process than exists in reality. Kinetic 
energy and pressure loss coefficients are the classical aero- 
dynamic parameters, but they do not accurately represent the 
performance of a cascade with injection. 

As a final observation or conclusion, it is noted that most 
studies correlate aerodynamic performance of film cooled blades 
with injection mass flow rate. This can be misleading to a 
designer since injection mass flow is dependent on the injection 
velocity and port area. It was found in this investigation 
that aerodynamic performance is more directly related to 
injection velocity than mass flow. The higher the injection 
velocity the larger the performance loss which is incurred. 

After the finish of the work on this grant, a discrepancy 
was found. It v?as determined that the turbine cascade configura- 
tion analyzed in this report was not the same as that of 
reference [16] . The stagger angle of the cascade used here was 
less. If the correct stagger angle v/ere used it would be found 
that the calculated aerodynamic performance curves for injections 
in the diffusion regions of the suction surface would be lower by 
approximately 0.003 than those shown in figures 19 and 20. 
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LIST OF SYMBOLS 


English Symbols 


A 


area 




a+bi 


B ,B 

s' V 


C 

c 


source-sink and vortex singularity influence 
coefficients for tangential or real velocity 
component at a point 

a general complex number or location in the complex 
plane 

source-sink and vortex singularity influence 
coefficients for normal or imaginary velocity 
component at a point 

complex potential function 

location of singularity distribution 


Cjj discharge coefficient for injection row 

C, combined defining matrix for use in relating the 

source-sink distribution to vortex distribution 



c 

P 


D 

e 

F 


f (X) 
f (z) 
G 

g 


coefficient of pressure 

specific heat at constant pressure 

end points of the line segment over which the 
singularities are distributed 

coefficient matrix of the minimization equations 

kinetic energy loss coefficient 

specified boundary condition function along 
the cascade body 

general real function 

general complex function 

singularity strength density function 

vortex singularity strength distribution 

total enthalpy 
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distance along the imaginary reference axis to 
the intersection point of the imaginary axis and 
an extension of singularity line segment 

the imaginary number, 

combined concentrated singularity strength of the 
cascade body as seen from the far stream boundary 
conditions 

combined singularity strength distribution, (s+ig) 

length of the sides of the polyhedron used to 
approximate the cascade body 

Mach niimber 

critical velocity ratio 

nximber of sides in the approximating polyhedron 
unit outward normal vector to cascade body 
static pressure 
total pressure 

static pressure at the outer boundary of the 
mixing layer 

constrained singularity strength density function 
for minimization 

dynamic pressure 

'radial vector from cascade 

gas constant for ideal gas equation 

radius or distance from a point in any direction 

cascade spacing 

source-sink singularity strength distribution 
static temperature 
total temperature 

streamwise coordinate along the cascade body, 
also trailing edge thickness 



V 


V 

w 

w 

X 


x+iy 

y 


z 


general velocity 

disturbance velocity due to the cascade 

weight or mass flow rate, also complex velocity 
due to a surface singularity 

complex velocity due to a singularity distribution 

real coordinate axis of the reference coordinate 
system 

general location in the reference complex plane 

imaginary coordinate axis of the reference 
coordinate system 

general complex number or location in the 
reference complex plane 


■Greek Symbols 


a 


B 


r 



A6 

n 

’^PA 


flow angle, also orientation angle of singularity 
line element 

cascade stagger angle 

onset flow contribution to the relation between 
source-sink distribution and vortex distribution 

concentrated vortex strength of the cascade body as 
seen from the far stream boundary condition 

circulation about cascade body 

ratio of specific heats 

finite incremental change in a variable 

boundary layer displacement thickness 

areas of influence of the source-sink and vortex 
singularities 

turning angle of the cascade flow 
orientation angle of singularity line element 
primary air efficiency 
thermodynamic efficiency 


54 



0 injection flow, angle to cascade body surface, 

orientation angle of the singularity line 
element and boundary layer momentum thickness 

X compressible flow correction factor 

Lagrangian multipliers of the constraints in the 
minimization equations 

? ratio of coolant flow to primary or mixing 

layer flow, also coordinate line along the 
singularity line element 

p static density 

total density 

Z concentrated source-sink strength of the cascade 

body as seen from the far stream boundary conditions 

0 cascade solidity 

<j) velocity potential 

i{/ stream function 

0 ) ■ pressure loss coefficient 

Other Symbols 


V 

(") 

Re( ) 

Im( ) 

^ hxk 
Subscripts 


gradient vector operator 
vector quantity 

real portion of complex variable 
imaginary portion of complex variable 
matrix with Z- rows and k columns 


BODY evaluated along the cascade body surface 

c coolant or injection flow quantity 

comp compressible flow corrected variable 

cr gas dynamic critical condition 

E,2 far downstream station or exit flow condition 

I far upstream station or inlet flow condition 
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matrix and equation indices 


inv 

inc 


m 


N,T 


onset 


sep 


t 


TE 

u 


x,y 


z 


CO 


inviscid flow variables 

incompressible flow variable, also incidence 
angle 

mixed flow condition downstream of an injection 
row in the mixing layer 

normal and tangential components of velocity and 
forces at points on the cascade body surface or 
anywhere in the flow field 

variables which are defined as constant throughout 
the flow field 

quantities which are associated with the injection 
induced separation model 

physical state variables which are evaluated at 
total conditions, also components of variables 
evaluated along the body surface 

quantities associated with the trailing edge of 
the cascade body 

flow conditions in the upstream or uhmixed region 
of the mixing layer before flow has crossed the 
injection row 

components of variables which are aligned v/ith 
the reference coordinate system 

quantity evaluated at any general point in the 
flow field 

quantities evaluated at stations far removed 
from the cascade row 
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FIGURED ONE DIMENSIONAL MIXING MODEL 
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UNIFORM FLOW 


FIGURE 4. CASCADE LOSSES WITH 
INJECTION 



FIGURE 5, INVISCID CASCADE PROBLEM 
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C,+iS 



FIGURES. SINGULARITY DISTRIBUTION 
' ■■ GEOMETRY 
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POTENTIAL FLOW 
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FIGURE 14; COMPRESSIBLE SOLUTION - TURBINE BLADE 
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FIGURE 17; PRESSURE DISTRIBUTION WITH AND WITHOUT INJECTION 
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At7p/-^p^^- CHANGE IN PRIMARY AIR EFFICIENCY 
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FIGURE 2L. MULTIPLE INJECTION ANALYSIS 
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FIGURE 22. DESIGN TURBINE. BLADE COOLING CONFIGURATION 



')7^^ANDi7p^ THERMODYNAMIC AND PRIMARY AtR EFFICIENCY 



FOR 
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AND?7p^THERM0DYNA 



FIGURE 24. AERO-PERFORMANGE OF DESIGN CASCADE 
FOR VARIATIONS IN INJECTION ANGLE 
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VARIATIONS IN INJECTION TOTAL 
TEMPERATURE 
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FIGURE 26; DESIGN CASCADE PRESSURE DISTRIBUTION 



APPENDIX 


MATRIX COEFFICIENTS 
' Singularity Influence Matrices 

The. problem is to calculate the velocity contribution that 
a distributed singularity has on a point. The singularity and 
the point have coordinate systems with varying orientations to 
a reference coordinate system. Resulting velocities are to be 
described in the point coordinate system. The general problem 
is shown in figure A-1. The reference coordinate system is x,y 
with points z, c^ and C 2 locating the point coordinate system 
(p,q) and the singularity coordinate The point and 

singularity coordinate systems are at an angle 9 and 4> with re- 
spect to the reference. All coordinate systems are complex. 

The prime notation indicates the coordinates of a point in the 
singularity or point of interest coordinate system. The surface 
singularity is distributed along the line segment of length £ 
which extends to' + » into and out of the plane. For a cascade, 
identical surface singularities occur at intervals of the cascade 
spacing, S, along a line parallel to the y axis. In this dis- 
cussion the simimed singularity of equation (31) in the analytical 
solution section is integrated to + <» in the out of plane direction 
and will be considered as a single singularity or velocity potential 
operating over the surface from c^, C 2 ; into and out of the plane. 


w = 



2 00th 


(z ’ -n) 
^ s 


(Al) 


where 


k{n) = s(n) + igCn) (^2) 

Unit vectors aligned with the real axis of the point and singu- 
larity coordinate systems may be written in complex notation. 
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. point unit vector = 
singularity unit vector = 
Constant Strength' Singularity ; 
larity strength is written as. 


(A3) 

(A4) 

For this case, the singu~ 


k(c) = k(ri) = constant 

The velocity effect in the reference coordinates is found as 
before in the analytical solution section 


^2 

w = 2 I coth dc 

^1 

sinh{TT (z-Cp)/S) 

w - - 2fc 5^* nn > 


w (z") 


v - iv = e^® w(z) 

P q 


w(z”) 


2k e^^e^*^ Jin 


sinh{7r {z-.n 2 )/S) 
^sinh(ir (z-n^)/S) ^ 


(A5) 


(A6) 


(A7) 


Equation (A7) applies for the general problem where the point of 
interest is removed from the singularity, but when the point lies 
on the singularity surface a Cauchy integral technique must be 
used to properly evaluate equations. When 


z’ (n,?) = z"CP/g) = 0+Oi 



then equation (Al) becomes 


w(z) = 2 ^ coth (“ “) 

Following the same procedure as before the integral form is 
written as 


+ £./2 

w = 2 ™ e^*^ / coth (- ~) dn 

-%/2 

The integration is carried out piecewise to isolate the singu- 
larity cross-over point at n = 0. 

“6 +6 

w = 2 e^® [/ coth(::^)dn + / coth (‘:^) dll 

-il/2 -5 


+ / coth(-~)dn| 

where +6 -define a small region around the origin. 

Since the hyperbolic cotangent is an odd function along the 
real axis^ and the distance —6 to —a/2 and +6 to a/2 are the same 
the first and third integrals will exactly cancel. 

+ 6 

w = 2 ~ e^*^ / coth(-|^)dn 

-5 

In the neighborhood of the origin the coth (z) may be replaced 
by 1/z and by letting u = -irn/S, the integral is simplified. 

+ 6 

w = - 2k e^*^ ^ ^ 

-6 
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This integral can be evaluated by using an indented contour to 
arrive at the Cauchy Principle Value 

— p +p +6 

w =- - 2k ■ f i du + / I du + / I du j 

-S ~p -p 

The first and third integrals cancel leaving 

w = - 2k {~iri) 

w(z”) = w(z) 

w(z") = 2-nrki (A8) 

The influence matrix is formed by using the above formulas 
(A7-A8) to calculate the summed effects of all the body surface 
singularities on each of the central points on the body. The 
final form of the constant singularity strength influence matrix, 
A, is a n X n matrix with n being the nvimber of singularities and 
control points. A single element of the matrix can be written as 

n 

aij = f 2iri - I 2 e^®i 

S-1 

a+i 


sinh(iT (z^;-C2. )/S ■ 

^sinh (tt (z^-c^. )/S^ 


and is the coefficient of the strength k^. 

In the matrix form of equation (54) and (55) of the analytic 
solution section for calculating velocities and the source-sink 
strength distribution (Sj^> , the source- sink influence matrices 
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Re [A] 


(AlO) 


ij nxn 


nxn 


xj nxn 


-IM [A] 

nxn 


(All) 


Linear Strength Singularity ; The. development of this coeffi- 
cient matrix follows closely the development of the constant 
strength matrix. The only change being th e form of strength 
function k in equation (A2) . 




k(n,l- k(n, ) 


k(n) 

= k(nj_) 

n - -n (n - n-, ) 

^ 2 ^ X 



= k(nj^) 

4- ACn ~ Tij^) 

(Al2a) 

kCn) 

= K(c^) 

+ A (z-c^) 

(A12b) 


Integrating the single surface singularity to find velocity 
contribution at a point gives 


w(z) 


2 Ih(n 3 _) / coth dTi 

^1 

t A / (n— n^) coth dn] 


(A13) 


Let 


D = / 


coth(:L4:^) dn 


„ sinh[ir {z-c^)/S3 

_ ^ r ^ 1 

7T sinh Eii (z-c^)/S3 ^ 


(A14) 


then 

w(z) = -2e Lc^] D - |2- / n coth dn > 
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(A15) 



Direct evaluation of the remaining integral in equation 
(A15) , results in an infinite series solution in which the coef- 
ficients of the terms are Bernoulli numbers. Rather than use this - 
complicated solution form, an integral is formed in a piecewise 
manner using the defining complex series for the hyperbolic 
cotangent function 

I = I {/ n coth ) dn) ■ 

n-. 


Expanding in series form 


I g { / n (tt z.-n) s ^ ^ 2 . , — 

T1-I IT +1T (z'-n) /S 


4Tr^+ir^(z'-n)^/S 9Tr^+ir^ (z ’ -n ) 


+ ...]} 


I 


^2 

/ 

^1 



« ^^2 

dn + J / 

a=l nj_ 


^ a„ 

(z'-n)+(jS)^ 


(A16) 


Equation (A16) can now be integrated terra by term to give a 
recursive relation for the integral 


Z-Cp » 

I = (c^-c,) - z 5,n ( " " c ' + I' 

1 a=l 

+ 2jS[tan - tan ] 

(jS)^+(Z-Cp)^ 

- z £n ( 5 ^)} (A17) 

(jS)"=+(z-c 


The index j corresponds to the number of bodies in the cascade ' 
row on either side of the point of interest being considered in 
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the velocity calculation. As the quantity, jS, becomes larger, 
the series terms contribute less to the velocity. In this solution 
the series evaluation was set up to include terms for the con- 
tributions of all cascade blades within distance of ten blade 
chords . 

Combining all the terms in equations (A13) , (A14) and (A17) , 
the disturbance velocity becomes 


w{z'') = -2e^^ e ^"^[{kCc^) - 


[k(c2)-k(c^) ] 


c^} D 


[k (C2)~k(c^) ] 




I] 


Gathering terms on the values of k 

^(gtr) _ k(c^) {-2e^® e [D + — -y (c^D +1)]} 

2 1 

c, D+I 

+ k(c2) {2 e^® ■■ ^ -3 } (A18) 

2 1 

Again the effect of a surface singularity upon its own control 
point must be considered. The development is straightforward. 

It uses the evaluation of D from the constant strength analysis , 
and direct substitution into equation (A17) to find the value 
of I. - 


at z’ = 0: D = — iri 

C3 

at z’=0: I = -Jl+ ^ 2 [-il+2 j Stan ^ (-g^) ] 

j=l ^ 

Stibstituting into equation (AlS) with c^ = -1/2 and C 2 = +H/2 
gives 


wCs'*) = k(-£/2) {rri+[2-4 I (-1 + 

j=l 


2jS , -1 , S, . . , 

tr '2js> > ' 


CO 

t k(il/2) {-Tri+E-2+4 I (-1 + ^ tan“^ (^) ) ] } 


(A19) 


j=l 
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To formulate the linear matrix B, the coefficients for the 
singularity strengths at the end points of the approximating 
flat elements must be determined. The linear distribution 
required that the- first and ia'st elements have the same 
strength (i.e, • This causes several special 

cases to be added to the coefficient equations . For j 7 ^ m^ 
and j 7 ^ 1 , the coefficient of k^ is 


^ jm 


ie -i<(). c.D.+I. 

e ^ {D t 

3+1 3 


ie 


+ 2 e 


■“iifc • T c. ,D. , +I . 

m ’' 3-1 

C .-C . ^ 

J D-l 


(A20) 


For- j 7^ m and j = 1 

ie^ -ii^T c,D,+I, 

Im ^1 (C2-C1) 

ie ~i(j) c D +I 

^n^nn n^ 

+ 2 e e {- 

^n+1 ^1 


(A21) 


For j = m and m 7 ^ 1 


mm 


j=l 


m 


oo * 

= {Tri + [2-4 I (-1 + tan’"^(^))]} 


le -Idi T C ,D ,+I T 
. m „ ^m-lr m-1 m-1 m-1-, 

+ 2 e e 1 } 


c -c , 
m m -1 


(A22) 


For j = m+1 and m ^ n. 


ie 


m +1 m 


= ~ 2 e 


m 


-icj. 


m +1 


c , , D , +I , T 
f-n J. n'+l m+ 1 . 


+ {—Tri + [— 2 + 4 ^ 

j=l 


(-1 + 


2 jS 

m 


-1 


m 


tan (A23) 
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For j = m+1 and m = n 


In 


-2e e {D^ + 


+ {-iri + [-2 + 4 ^(-1 + 1^. tan ^ ) U 


-l/n 


(A24) 


For j = m and m = 1 


1 = {iri + [2“ 4 ^ (-1 + 

j=l 




ie^ ~i(j) c D_+I 
+ 2e 1 e ^ 

°n+l" n 


(A25) 


With equations (A20) thirough (A25) , the coefficient matrix B 
may be formed for use in matrix form of (54) and (55) of the 
analytical solution section for the calculation of body velocities 
and determination of the vortex strength linear distribution (^j^) - 


[A' ] = ReliB] 

^ij nxn nxn 

IB ] = -ImliB] 

^ij" nxn nxn 

Minimization Equations Matrix Form 

Recalling the general form of the equation system for finding 
the minimized singularity density, the formula needed for 
developing the details of the coefficient matrix may be written 
down. 


(A26) 

(A27) 
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Q *■ Q t 9^2 * • • • 


93 t / ^ 2 ^ 


{[A_ ] 

H[3^] 

+ EC.j^] 

t^k^ . . 

“i Ixn, 

nx-1 

nxn 

nxl 

[A ] 

V. 

X Ixn 

[gvl J + 

^ nxl 

Xi{[B ] 

TEi 

E3.] 
Ixn : 

[B ] 

TEi Ixn nxn 

I5k3 ^ 

nxl 

f^v 

TEj 


nxl 


nxl 


" '^onset ^ 

" = 0 

Taking the derivative of Q with respect to each of its independent 
variables and setting the result equal to zero, generates the 
minimization matrix equation. The derivative of Q with respect 
to g^ is in expanded form, 

IgT ^^s-, ^Ij ^^l‘^‘^ll^l'*’^12^2 + ••• ^Ij^j ^ •••‘^In-^n^ 

3 -*- 


+ ... 2Ag^ Cjj (3j+Cjj^g^+Cj252 + ••• ^jn^n^ 


33^3 


+ ... 2Ag^ C^j( n'*‘‘^nl^l‘*‘^n2^2’ ''' *** *^nj^j ■*■*** ^nn^n^ 


+ 2A • g. + C . + B ^ . 


n3"D 

+ ... B 


s j j 

TEj .. 


+ ...B C.+B )+\-(C.-Ct.) = 0 (A28) 

®TEn “3 Vjj,. 2 nq iq 


The derivative with respect to is 
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®TE1 TE2 


6 +...B 

2 ®TEj “ ■ 


^ ^ 12^2 “^ * • • 

■^ ••• *^j2^2 ■^ ••• 


+ . . . B (C T g. + C + . . . C g ) 
Sm„ nl ^1 n2^2 - nn^n 

xi^n 


5 + B g_ + . . . B g_ 

'^TEl ^ ''tE 2 ^ '^TEn 


- '^onset -“> = ° 


The derivative with respect to ^2 is 


“ ^^nl ' ^11^ ^1 ^^n2 " ^12^ ^2 




- (e^ ~ =0 (A30) 

After collecting multipliers of the independent variables and 
taking constant terms to the right hand side of the equations , 
the general coefficients are written. The coefficient of g^^ in 
equation (A28) when j 1 is 


2 (A C, . Ct . + A C„ . Ct + . . . A . C. . C . . 
' 1] lx S2 23 2i S3 X3 33 

+ ... A C. C’) 

®n ^ 

and when 3 = i , coefficient of g^^ becomes 


• 2 (A C?. + A C?. + . . . A . . + . . . A C^.) 

I3 S2 23 S3 33 n3' 
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The coefficient of in equation (A28) is 


(B Ct . + B . + . . . C . + B , ) 

I 3 ®TE2 ®TEn ^TEj 


TEl 


and the coefficient of 


- ‘^ij> 


Finally, the equation constant is 




• + B„) 

s^ n 3 n 


For equation CA29) the coefficient of is 


(B . + B . + ...+ B C . + B ) 

®TE1 ®TE2 ®TEn ^TEi 


The coefficients of X^ and are 


( 0 ) 

and the equation constant is 


^’^^®TE 1 ^®TE 2 ’** ^n^ ^onset 


For equation (A30) , the coefficient of is 


«=ni - =li» 


The coefficients of X^ and are 


( 0 ) 


and the equation constant is 


(Bi - e„) 


sin(e^g-a) ] 
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Using the coefficients and constants, a. system of linear 
equations results which is solved for the minimum vortex strength 
density. 


3Q 

3g, 


coefficients 


3Q ~ 
9go 


coefficients 


3Q 

3g. 


~ coefficients 


n 


'n 


SQ 

ax. 


coefficients 


10 . 

ax. 


~ coefficients 



constant 

constant 


constant 

constant 

constant 


(A31) 
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FIGURE A-L SINGULARITY INFLUENCE PROBLEM 
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